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ABSTRACT 

Here we revisit line identifications of type I supernovae and highlight trace amounts of unburned 
hydrogen as an important free parameter for the composition of the progenitor. Most 1-dimensional 
stripped-envelope models of supernovae indicate that observed features near 6000 — 6400 A in type I 
spectra are due to more than Si II A6355. However, while an interpretation of conspicuous Si II A6355 
can approximate 6150 A absorption features for all type la supernovae during the first month of free 
expansion, similar identihcations applied to 6250 A features of type Ib and Ic supernovae have not 
been as successful. When the corresponding synthetic spectra are compared to high quality time-series 
observations, the computed spectra are frequently too blue in wavelength. Some improvement can be 
achieved with Fe II lines that contribute red-ward of 6150 A, however the computed spectra either 
remain too blue, or the spectrum only reaches fair agreement when the rise-time to peak brightness 
of the model conflicts with observations by a factor of two. This degree of disagreement brings into 
question the proposed explosion scenario. Similarly, a detection of strong Si II A6355 in the spectra of 
broad-lined Ic and super-luminous events of type I/R is less convincing despite numerous model spectra 
used to show otherwise. Alternatively, we suggest 6000 — 6400 A features are possibly influenced by 
either trace amounts of hydrogen, or blue-shifted absorption and emission in Ha, the latter being an 
effect which is frequently observed in the spectra of hydrogen-rich, type H supernovae. 

Keywords: supernovae: general 


1. INTRODUCTION 

Studying the evolution of a supernova’s spectral en¬ 
ergy distribution is a formidable task on account of ex¬ 
cessive line blending that, for all intents and purposes, 
erases most traces of a reference continuum level. Specific 
line transitions that contribute to conspicuous features in 
SN la spectra have been mostly identified, whereas the 
identification of weak spectral features for most SN I su¬ 
pernovae remain uncertain (see Table I). Line blending 
on the scale generated by the large expansion velocities 
(> 10"* km s“^) thus creates a number of hurdles for 
supernova spectroscopists who are interested in extract- 
ing information about the comp osition of the ejecta (cf. 
Gaposchkin||I936 Popperp937 ). 

In terms of a taxonomy, or a observed spectral se¬ 
quence, supernovae of various luminosity classes are di¬ 
vided into two central categories, types I and H, depend¬ 
ing on whether or not th e spectrum contain s conspicuous 
signatures of hydrogen (Minkowski 19411. For type H 
supernovae (SN H), signatures of hydrogen Balmer lines 
come in the form of distinct P Cygni profiles. Depend¬ 
ing on the decline in the light curve, type H super¬ 
novae can be further sub-classified as SN HL for a linear 
decline and SN HP for sustained brightness before an 
eventual decline. However, drawing distinctions between 
SN HL and HP is less obvious co nsidering an existing 
continuum of observed properties (Sanders et al. |2014 


.|2014[|Peicha fc Prieto|2U15||Call et al. [21)15 

ti et aL| |2015p . 
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Broadly speaking, hydrogen-poor SN I supernovae can 
be separated into three subclasses, SN la, Ib, and Ic. 
In addition, the spectra of SN Hb evolve from SN H 
in appearance, with weak Ha, but later transition to a 


spectrum w ith lines of helium similar to SN Ib (Filip- 


penko 1988). Due to similarities between some caught- 


early SIN ib to SN Hb at later times, some SN Ib may 


be ca ught-l ate SN Hb (|Folatelli et al.||2014[ see also Liu 
et al.|[20T5 |. 


Types Ib and Ic are generally distinguished by the pres¬ 
ence o r abse n ce of obvious he l ium sig natures (Wheeler] 
et al. 1995 Dessart et al. 2015bI. Additionally, 


thermonuclear SN la reveal conspicuous signatures of 
Si II A6355, while the latter varieties of SN I core-collapse 
supernovae do not (|Filippenko| |1997|); i.e. stripped- 


envelo pe subclasses of SJN Ib/Ic are dehned by Filippenko 
(19971 and previous authors as not having a conspicuous 
doublet signa ture of Si II A6355, w hich is a line detected 
for all SN la ( Wheeler et al.|[l9% ). 

For a comparison of these types of supernovae, in 
Figure we plot the SN Ib 1984L, SN Hpec 1987A, 
SN Ic 1987M, and SN la 201 He and note the similar¬ 
ity of 6250 A features between SN H-peculiar 1987A 
and SN I c 1987M is strikin g. However, in the similar 
Fig. 1 of Filippenko (1997), weak absorption features 
near 6250 A are labeled as Si II for both “SN Ib” 1984L 
and “SN Ic” I987M, implying the spectra harbor a lone 
and weak feature of silicon when this feature may in¬ 
stead be in fluenced by trace amo unts of hydrogen for 
some SN Ib ( James fc BaronpOlO ). 

A precise d efinition of SIN Ib/Ic has since vari ed over 
the years (see Fryer|20(M{ Gray &: Corbally|2009 ). Some 
interpretations assume hydroge n is absent tor SN I “by 
definition of being type I” (cf. Maguire et al. 2010 El-| 
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Type I and Type II Supemovae 



5400 5800 6200 6600 

Rest Wavelength [A] 


Figure 1. Comparisons of spectra from SN 1987A soon after the explosion on UT February 24 , 1987 JPun et alT]|1995ll: S N 1984L at 
approximately day +8 relative to S-band maximum light l|Wheeler &: Levreauit||1985||Doggett & Brancli|iy85||T'svetkov|i987|l; SN 1987M 
on UT September 28. 1987 and UT Oct ober 2, 1987, approximately days + 9 and -4-iij. respecti vely, relative to .tf-bancl maximum light 
< |Filippenko et al.|1990|[jeffery et al.|1991[ |; and SN 201 Ife near maximum light l |Pereira et al.|2013| l. Labels for SN 1987A and 201 Ife denote 
ions for which detections have been verified, while labels for the 6250 A features of SN 1984L and 1987M represent unproven candidate 
ions. For SN 2011fe, 1984L, 1987A a single epoch is shown, whereas for SN 1987M two epochs are shown (+9, +13). 
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dridge et al. 2013), while lone detections of Si II A6355 are 
deemed robust tor SN Ib/Ic, and therefore accurate for 
measurement (s ee, e.g., Taubenberger et al. 2006 |Mod- 


jaz et al. 2014). Line identifications for 6000—6400 A 


features in the spectra of broad-lined SN Ic (BL-Ic) and 
type I/R super-luminous supernovae (SLSN) have also 
been debated. Most notably, estimates of Si II and Fe II 
conflict when 6000 — 6400 A spectral features are in¬ 
terpreted to be dominated by Si II, while model spectra 
producing mostly signatures of Si II are consistently too 
blue in wavelength when compared to observations (cf. 
Mazzali et al.|[2000[ Inserra et al. 2013[ Lyman et al. 


2044 ^ 


The case for detectable hydrogen i n the SN Ic 19941 


was d i scusse d by Wheeler et al. (1994), while Clocchiatti 


et al. (1996) claimed detections of weak lines ot He I as 


well. I ^th i nterpretations were later favored by [Parreiit 


et al. (2007) who used SYNOW to show that a mix of both 


photospheric and higher velocity H I A65 63 offer better 
altern atives to the unm a tched models of IBranch et al.l 
(2006) and Sauer et al. (2006), both of which produce 
strong signatures of Si II. The feature nearest to 6250 A 
has also b een proposed to be a c omposite of C II, Ne I, 
and Si II (Hachinger et al. 2012). However, the abun¬ 
dance of Ne I needed to match observations is reportedly 
too high to be physically consistent with the final abun¬ 


dance s of candidate progenitor systems (Ketchum et al. 


2008). 


Dessart et al. (2015b) recently produced non-LTE spec- 
trum synthesis calculations, stemming from “hydrogen- 
poor” (Mjj ~ 10“^ — 10“^ Mq) and “hydrogen- 
depleted” models (M// ^ 10“® M©), and have shown 
that a Si II/Fe II-weighted prescription for 6250 A fea¬ 
tures of SN Ib/Ic may provide a more plausible alterna¬ 
tive to treating weak signatures of Ha as a free parame¬ 
ter. These models also indicate that if the progenitor star 
has above 10“^ Mq of hydrogen within a thin shell on 
its surface (i.e. not physically detached), then the cor¬ 
responding post-explosion Ha profile near 6250 A will 
resemble that of SN Hb so long as the surface abundance 
of hydrogen is high enough. However, the rise-time of 
these SN Ib and SN Ic models are roughly twice as lon g 


as what is typically observed (cf. Wheeler et al. 2015), 
which questions the models ability to represent ^'stan¬ 
dard compositions”. 

Finally, recent observations of the SN Ib 20I4C reveal 
evidence for interaction with a hydrogen-rich circumstel- 


to explosion dMilisavljevic et al. 20151 see also 

Chugai 

& Chevaliei'l 2006| |Yan et al. 2015 

Moriya et al. 

2015). 

bhould we then expect trace amounts ot hydrogen to re- 


main within the outermost layers of the progenitor for 
some of these type I core-collapse events? 

In this work we examine the latest model spectra for 
type I supernovae and confront them with observations. 
In §2 we highlight some current pitfalls when invoking 
detections and non-detections of select atomic species. 
In §3 we overview the ability of various explosion mod¬ 
els to reproduce 6000—6400 A features in the spectra of 
SN la, Ib, Ic, and more luminous events. In §4 we sum¬ 
marize our findings. 

2. SPECTRAL LINE IDENTIFICATION 


Line blending is so severe for supernova spectra that 
it permits broad consistency between a synthesized spec¬ 
trum and a majority of observed spectral features. These 
include absorption and emission line signatures that can 
be either conspicuous or weak below limiting noise. Some 
of these features can be approximated to a few lines, 
while others sustain a complex blend from first light to 
well after a year since explosion. 

Whether testing a benchmark model, or tracing ob¬ 
served features with a prescription of candidate lines, the 
interpretation of the observed spectral features, and the 
related empirical measurements, depend on the quality 
of the full fit from the supporting model. Moreover, if 
a model spectrum produces all but a single (or few) ob¬ 
served feature(s), the model is judged to be significantly 
incomplete. The same is true if the model produces spec¬ 
tral signatures that are not observed at the same epoch of 
the observation. Considering the relationships between 
the composition of the ejecta and the energy throughput, 
when a model’s light curves are in good agreement with 
observations, the model spectra must also reflect obser¬ 
vations before a model is sufficiently matched, and vice 
versa. 

Since boundaries between spectral line signatures are 
not well-defined for most supernovae, empirical false- 
positives are readily available within an underspecified 
parameter space of “detections”. Consequently, detec¬ 
tions and non-detections of various ions and spectral lines 
do not follow classical definitions. For example, a non¬ 
detection of the strongest line from an ion at 5500 A 
does not necessarily imply a lack of significant spectral 
contamination from either this or other ions at 5500 A 
and other wavelengths. This can be made worse by low 
signal-to-noise ratio data and inadequate spectroscopic 
follow-up. 

A common example is singly ionized silicon. This ion 
is considered detected in SN la optical spectra when a 
number of spectral lines are conspicuous above noise and 
additional blending from adjacent lines. However, addi¬ 
tional signatures of Si II AA3858, 5972 are not definitively 
detected for SN Ib, Ic, and SLSN-I. 

See Figure where we hav e plotted photosphe ric phase 
spectra of th e SN la 201 Ife (Pereira et al.||20I3 ) and the 
SN Ic I994I (|Filippenko et ar|I995) m terms of Si II and 
H I line velocities CI(v/c). For SN 20IIfe, its spectra 
are reasonably influenced by so-called conspicuous lines 
of Si II. In comparison to SN Ib/Ic, a lack of obvious 
Si II A5972 is to be expected since the stronger Si II A6355 
would already be a weak signature compared to that for 
SN la; i.e. a purported identification of Si II A6355 so 
far has only one wavelength region by which to constrain 
the associated model. 

Other ions with few strong lines, such as H I, C II, 
C III, Na I, and Si IV, face similar obstacles during pho¬ 
tospheric phases as well, while a signature of Na I D 
is arguably absent for SN la, minimal or absent for SN 
Ib, and significantly present for SN Ic (Dessart et al. 


2012 2014b 2015b). Hence, studies of ejecta composi¬ 
tions benefit when tne number of candidate lines, is either 
effectively maximized or close to the number of lines that 
are suspected to be present ( Baron et al. 1996 Friesen 
et al.|[20l4 Vacca et al.||20I5 ). 


These are important points to consider when interpret- 
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H-Balmer Frame Si II Frame 

Ha Hp U6355 5972 5051 4130 



- 30-20 -10 0 -20 -10 0 

Projected Doppler Velocity [km s'^ 1000'^] 

Figure 2. Plotted above are projected Doppler velocities 0(v/c) in the frame of the strongest H-Balmer lines (left) and Si II lines fright). 
The light g rey bands ar e meant to indicate the overlap bet ween common signa tures of a given ion. Spectrum references: SN 19941 ||Modjaz| 
|et al.| 20T^; SN 1996cb Matheson et al.|2001[ |; SN 2011fe ( [Pereira et al.|201^ . Model references: (D14) PDDEL4, day —6.6 for SN 2Ullte 


by juessart et ah] | |2014a 





























5 


ing type I spectra given that freshly synthesized ejecta 
are potentially contaminated by unburned and distin¬ 
guishable traces of H, He, C, O, Si, Ca, and Fe-rich 
progenitor material near and above the primary volume 


of line formation (a so-called “ph otosphere,” see Hatano 


et al. 1999 Branch et al. 20051. Additionally, an ob¬ 


served spectral signature associated with a particular ion 
at photospheric velocities can be influenced by an iden¬ 
tical ion, in addition to other atomic species at si milar 
or relati vely higher velocities (Marion et al. 2013) (see 
Fig 3. of Blondin et al.|[ ' 


For example, Sasdelli et al.|( 2014) claim non-detections 
for all signatures of carbon in the early spectra of the 
peculiar SN la 1991T without providing an explanation 
for a tentative detection of C III A4649. This identifi¬ 
cation is also consistent with kinetic energies deduced 
from detectable lines of Fe 111 and a tentative detection 
of Si 111. For the computed spectrum at 4500 A, a fea¬ 
ture is missing exactly where it is present for SN 1991T, 
1997br, and 1999aa during early pre-maximum epochs. 
Therefore, abundance determinations of unburned car¬ 
bon would benefit from knowing the subset of species 
giving rise to weak 4500 A f eatures in the spectra of 


SN 1991T/1999aa-like events (Hatano et al. 2002 Par- 


rent et al^[2011[ 2014|). 

Detections of select atomic species are most appar¬ 
ent when line-blanketing effects are strong. A clas¬ 
sic example is the absorption trough of Ti 11 observed 
for SN 1991bg and 1999by on account o f overall lower 
ejecta temperatures (Doull & Baron||2011 and references 
therein). By contrast, less conspicuous detections of Ti II 
will be blended and relatively ambiguous against a back¬ 
drop of signatures from other iron-peak elements. This is 
the root of the problem when sub-classifying supernovae 
according to their s pectral features . 

As an example. White et al. (2015) utilize non¬ 
detections of the Ti II absorption trough observed blue- 
ward of 4450 A in an effort to disti nguish peculiar 


SN 2002es- likes with conspicuous Ti II (Ganeshalingam 
et al.|20i2 ) from SN 2002cx-likes found without the clas¬ 


sically seen trough of Ti II absorption (Li et al.|2003 ). It 
is important to note, however, that despite clairns of non- 


likes are not likely void of Ti II (see 

Foley et al. 

2009 

Kromer et al.||2013 

Stritzinger et al.||2 

014). 


3. MODEL SPECTRA OF TYPE I SUPERNOVAE 
NEAR 6000-6400 A 


In this section we examine SN I synthetic spectra found 
throughout the literature and assess respective goodness- 
of-fits when compared to observations. For each model 
spectrum shown, we have attempted to preserve the flux 
scaling for comparisons between observed and synthe¬ 
sized data in the original works, and we have normal¬ 
ized the spectra with respect to the 6000—6400 A pro¬ 
files where appropriate. (This underscores the need for 
releasing all published model spectra and photometry.) 


3.1. Thermonuclear SN la 

Plotted in Figure are observations of SN la sub- 
types compared to the computed spectra from var¬ 
ious spectrum synthesizers at wavelengths spanning 
5200 — 7000 A. All model spectra, generated from semi- 
empirical to full non-LTE radiation transport solvers, are 


Si II in Type la Supemovae 



Rest Wavelength [A] 


Figure 3. Comparisons between computed and observed s pectra 
of select SN la. Spectrum ref erences: SN 1994D (br own). 
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capable of providing close match to the feature centered 
about 6150 A through prescriptions based on Si II A6355. 

Specihcally, all absorption and emission components 
are well-centered with the data, while the curvature of 
the profile wings are relatively well-matched. This is also 
in spite of either a “non-standard” distribution of Si II, 
or two distinct components of photospheric and detached 
Si II A6355. That is, precise and accurate predictions are 
achievable for the red and blue wings, which are both 
sensitive to mean expansion velocities and the radial de¬ 
cline in density. 

To further appreciate the limiting accuracy of 6150 A 
features in SN la spectra as having been identified as 
Si II A6355, consider that an intrinsic difference of 500, 
1000, and 2000 km s“^ in terms of line velocities would 
be reflected by an underlying shift of ^ 10, 20, and 40 A, 
respectively. That is, apart from opacity effects, signif¬ 
icant differences of 2000 km s“^ can give rise to 40 A 
blue-shifts that are of similar order to redshifts seen for 
6150 A features during the first month of free expan¬ 
sion. Despite these differences, model spectra for SN la 
consistently reflect a suite of observations. 


When confronted by events such as t he SN la 2012fr 


(Childress et al. 2013 

Maund et al.|2013 

), where the con- 

spicuousness ot a hig 

ler velocity signature ot bi 11 At)355 


is sensitive to the cadence of follow-up observations, 
approximate consistency with composite features near 
6150 A is relatively large for SN la. In other words, 
the relative location of high velocity Si II would not be 
well-constrained without the accompanying time-series 
spectra. As we will discuss in §3.4, this window of con¬ 
sistency is even greater for broad-lined SN Ic (BL-Ic) and 
SLSN. 


One reason a model SN la might not match near the 
emission component of the Si II is when the model does 
not center C II A6580 at similar projected Doppler ve- 
locities as Si II during the first month after the e xplosion 
(Taubenberger et al. 2011 Pereira et al.|2013). Other 
lines that impact this region for SN la include S II A6305 
and lines of Fe II near maximum light and thereafter. De¬ 
spite these additional biases on direct measures of Si II 
features, the identification of Si II A6355 as a dominant 
line in the spectra of SN la is robust for a variety of 
normal and peculiar events. 


3.2. Core-collapse Supernovae 

Bas ed on the recent spect rum synthesis calculations 
from iDessart et al.| (|2015b), who utilize a modified 


version of CMFGEN rtHiilier & Dessart 


1 


20121, stripped- 


envelope models that produce weak signatures of 
Si II A6355 may avoid both consistent mismatch and the 
need for trace hydrogen by relying on con tribution from 
Fe II, as well as He I in some instances (Dessart et al. 


20121. In this case, Fe II effectively moves a signature 
dominated by Si II toward redder wavelengths, thereby 
making a blend of Si II and Fe II a strong contender for 
6 250 A features of SN Ib /Ic spectra. 


Dessart et al. (2015b) also draw association with the 
observed spectral sequence of SN lib, Ib, and Ic, but 
without comparing any of these models to observations. 
Notably , the SN Ib 6p5Axl model of [Dessart et al.| 
(2015b) has a rise-time of 40 days. This rise-time is 
nearly twice as long as those typically observed for SN Ib, 


while the SN Ic 5pllAxl from Dessart et al. (2015b) has 
a ri se-time that conflict s with most observations as well 
(cf. Wheeler et a l.|20T^ . The “Bmf5p 09” series of mod¬ 
els presented by Dessart et al. (2012) also exhibit rise- 
times that conflict with observations by as much as 30 
days. 

Since the modified version of CMFGEN used for super¬ 
novae has not been made publicly available, we are un¬ 
able to test areas of improvement f or any of the pr o posed 


models . The model spectra from Dessart et al. (2012 
2015b) have also not been made available tor comparative 


st udies, which was don e for the SN la models computed 
by Blondin et al.j (|2015|) . Therefore we take the liberty of 
digitizing model spectra of interest to plot them against 
observations in Figures]^ — 

3.2.1. “SN lb” 2007kj 

In an effort to distingu i sh SN Ilb, Ib, and Ic from 
their spectra, Liu et al. (2015) recently conducted a 
comparative study ot “the largest spectroscopic dataset” 
of stripped-envelope supe rnovae. Using the Supe rnova 
Identification tool, SNID (Blondin & Tonry 2007), and 
pseud o equivalen t meas urements ot several spectral fea¬ 
tures, [ Liu et ah (2015) classify SN 2007kj as a SN Ib. 


iVIodjaz et al. 1 201^ previously used SNID to classify SN 


2007kj as a SN Ib as well. 

However, comparison to the SN Hb 3p65A xI model 
spectrum pre^nted by Dessart et al. (2015b), plotted 
in our Figure would imply that SN 2007kj is instead 
either a caught-late or slow-evolving SN Hb. That is, 
in spite of having incorporated the large st set of SNID 
templ ates to classify SN H b, Ib, and Ic, [Modjaz et al.j 
(|2014|) and [Liu et ah (|20I5D are unable to type SN 2007kj 
as a SN Hb. 

Furthermore, if the SN Hb 3p65Axl model had been 
built for SN 2007kj, the match to SN 2007kj would sug¬ 
gest the progenitor had ^ 10“^ — 10“^ Mq of hydrogen 
on its surface at the time of explosion. It is also impor¬ 
tant to note that the 6250 A feature on day —2 for the 
SN Hb 3p65Axl model is not contaminated by Ha. This 
is in spite of a similar 6250 A feature on day —14 (and 
earlier) where an Ha signature is present for the model. 
By maximum light, this feature is pr edicted to be dom¬ 


inated by profiles of Si II and Fe II (Baron et al. 1995 
Dessart et al.||201fib ). 


3.2.2. SN Ib 1999dn 

The computed spectra of the SN Ib Bmf4p41x4 model 
are also associated with a p rogenitor having ^ I0~^ Mq 
of hydrogen near its surface ( [Dessart et al. 2012). See the 
bottom-left panel of our Figures and The difference 
to note between the day —2 spectrum of the aforemen¬ 
tioned SN Hb 3p65Axl and SN Ib Bmf4p4Ix4 models 
is that the shape of the 6250 A feature in the latter is 
not SN H-like, nor is its spectrum “Ha-positive”. This is 
in spite of the hydrogen present in the ejecta, while the 
6250 A feature is also blended with He I. 

Looking at the middle panel of Figure we find a 
promising match between the SN Ib I999dn at day —10 
to the SN Ib Bmf4p4Ix4 model at day —16. However, 
this model spectrum is not well-matched since the He I 
features near 5650 A and 6900 A are significantly too red 
in wavelength. Increasing overall kinetic energy for the 















































































Scaled Flux + Constant 


7 



Rest Wavelength [A] 


Figure 4. Comparisons between comp uted spectra, SN 20 07kj (|Modjaz et aL||2014||, and SN 1999dn (|Deng et aL|2000|l. The data have 
been scaled and normalized according to Jeffery et al ||| ^07ll. In tne tour corn er panels^ we plot time-series spectra ot tne SN Ilb 3p65Axl, 
SN Ib 6p5Axl, and SN Ic SpllAxl models pre sented b^Uessart et al.|||2015b||, and also include the earliest spectra of the “standard SN Ib” 
Bmf4p41x4 presented by |Dessart et al.| | |2oi2[ |. The rise-time, initial mass, and final hydrogen mass for each model are listed in each panel. 
Dotted lines denote when tne spectrum is contaminated by signatur es of Ho. Certain featu res are labeled by ions producing significant 
structure within the computed spectrum. JBIO PHOENIX calculation of| James & Baron| (|2010|l with and without hydrogen are shown as the 
dashed black and solid blue line, respectively. 


model, among other variables, can resolve this discrep¬ 
ancy. However, the issue that remains is whether or not 
this increase would result in a 6250 A feature that is too 
blue in wavelength. Even so, if a match is deemed secure 
for a model such as Bmf4p41x4, then this would imply 
the progenitor had ~ 10“^ Mq of remaining hydrogen in 
its outermost la yers, as a SN Ib. 


3.2.3. SN Ib 2008D 


Alternatively, James & Baron (2010) previously pre¬ 
sented PHOENIX calculations to show that unburned hy¬ 
drogen at high expansion velocities (in a shell) is rea¬ 
sonable for SN 1999dn. (This comparison is also shown 
in Figure ffl) In this paradigm, the detachment veloc¬ 
ity of the nydrogen shell is treated as a free parameter. 
With the shell displaced to 19,000 km s“^, the mate¬ 
rial needed in the atmosphere of SN Ib 1999dn to re¬ 
produce the 6250 A feature (M// < 10“^ Mq) is less 
th an 0.01 Mf7^< Mff < 0.2 Mq inferred for SN Bb 1993J 


by Baron et al. 


(19951. 


This mass estimate does not 
also conflict with the swept-up 5 x 10“^ Mq estimated 
from model comparison s to late-nebular spectra of SN Bb 
( Jerkstrand et al.||2015 |. 


Dessart et al. (2012|) compared their Bmf5p09x4 model 
corresponding to 10.3 days since explosion to the Febru¬ 
ary 2 spectrum of the SN Ib 2008D. In terms of the num¬ 
ber of days since maximum light, and with a rise-time 
of 36 days for the model, this corresponds to drawing a 
comparison between a model at day —26 and an observed 
event at day -1-4, a difference of one month. 

In Figure we compare the day —34 spectrum of 
Bmf5p09x4 (the earliest available) to the day —13, —15, 
and —17 observations of SN 2008D. Apart from the dis¬ 
cordant epochs, the model spectrum does not confidently 
resemble those of SN 2008D prior to maximum light. The 
data are also fairly noisy at these early epochs, however 
this particular model is not optimized for drawing con¬ 
clusions on the validity of Si II A6355 for SN 2008D. 

If both computed and observed spectra were those of 
SN la, the differences are enough to signify different 
subtypes that may or may not be from the same fam¬ 
ily of progenitors. Specifically, and compared to SN la 
where detections of Si II have been verified, the earliest 













































Scaled Flux + Constant 



Rest Wavelength [A] 


Figure 5. Same as Figure |4| except now for SN 2008D llSoderberg et al. 120081 IMod i az et al.|2009|l and SN 19941 (|Modjaz et al.|2014f. We 
have replaced the Bmf4p41x4 spectra with those of BmtopU9xl trom|uessafr~et al.| ||2U12||. 


computed spectrum of Bmf5p09x4 does not overlap the 
6250 A feature in the spectrum of SN 2008D near the 
date of explosion. One expects improvement through 
either stronger contribution from Fe II or lower mean 
expansion velocities attributed to Si II, or the model is 
flawed ab initio. 


3.2.4. “SN Ic” 19941 

For SN 19941, a so-called “standard SN Ic”, two op¬ 
posing interpretations prevail in the literature for two 
models that produce mis matched signatures of Si II com¬ 
pared to observations (cf. |Branch et al.|2006| |Sauer et al. 
z006|. See Figure|6|where we have over-plotted these two 
i ndependent results. 

Sauer et al.| (|2006[) predict the feature centered about 


6200 A in the photospheric phase spectra of SN 19941 is 
primarily a product of C II, Ne I, and Si II, with mostly 
Si II present near maximum light. Based on their series 
of model spectrm including the fit n ear maximum light 
shown in Figure 6 Sauer et al. (20061 claim a satisfactory 
match has been made. 


Branch et al. (2006) previously used SYNOW to show 
that a satisfactory solution can be obtained for SN 19941 
near maximum light when He I, O I, Ca II, Ti II, and 
Fe II are assumed to form within similar strata of ejecta. 


When Branch et al. 
SYNOW ht, the predic 


(2006) incorporate Si II into the 
ed ottset of the feature dominated 


by Si I I is n early identical to that later found by [Sauer] 


et al. (2006). Thus, while it is reasonable to expect a 
multncomponent prescription throughout the spectra of 
SN Ib/Ic, it is unreasonable to suspect a correct interpre¬ 
tation for over-prescribed C II, Ne I, and Si II when the 
model used to support this prescription is incapable of 


producing the complex evolution of the data (cf. Baron 
etay|19^ Ketchum et al.|[2008 1. 


Hach inger et al. ] 2012| ) later used the model of |Iwamoto 


et al. ijTUM ) to argue that Si II A6355 is a viable op- 


tion for the 6250 A feature in SN 19941, and that this 
signature of Si II can be contaminated by other species 
including Ha. However, despite offering to remedy the 
mismatch of model spectra with increased contribution 
from Hg, the hydrogen /helium-enriched C-|-0 models of 


Hachinger et al. (2012) r^ain dominated by ill-matched 
Si 11 A6355. See Figure and note the asterisk next to 
the orange H12* spectrum near 6200 A; this dip is the 
Hg coming through. 

I n the top of Figurej^ we have plotted the fit produced 
by Parrent et al. (2007|) with SYNOW where contribution 
from Hg was treated as a free parameter. It is important 
to reiterate that the assumptions of the SYNOW/SYNAPPS 
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Rest Wavelength [A] 

Figu re 6. Model spectra compared to SN Ic 19941 lIMo^jazetal. 
| 2014[l. Model references: (P07) SYNQW calculation offl-' arrent et al. 

with hydrogen treated as a free parameter; SOb ll sauer et all 
2UUfc)ll ' B06 lIBranch et al.ll2006ll: H12 models are frornlH^Einger 
et al.|(|2012|l' 4 he blue model, 1112, is the base SN Ic 941 model lb 
days alter explosion in the top of their Fig. 11 in green. The or¬ 
ange model, H12*, corresponds to the hydrogen-enriched SN Ib 
in the top of their Fig. 11 in red. The black asterisk is near 
the location of Ha coming through in this model. Like SOb and 
B06, both H12 models fail to match the minimum 6250 A in 
SN Ic 19941 as they too are reliant on significant contribution from 
Si II A6355. Vertical-dashed lines denote the rest wavelength of 
Si II A6355 (blue) and H I A6563 (red). 


model are not realistic. Therefore best guesses for a 
composition, without any prior information from more 
detailed investigations, cannot necessarily be considered 
superior on account of a well-matched output spectrum. 
However, the primary purpose of SYNOW-like tools is to 
quickly explore where detailed models are consistently 
unsuccessful. 

To compare wi t h the m ore recent “SN Ic” calculations 
of|Dessart et al.| (|2015b|), in Figure we have included 
tiine-series observations of SN 1994! and over-plot the 
earliest computed spectrum of SpllAxl at day —31 with 
the day +23 spectrum of SN 19941. Overall, one can see 
that the match is promising. However, the early evolu¬ 
tion of the 6250 A feature observed for SN 19941 remains 
unaccounted for by the model. The 6250 A feature at a 
particular epoch may be dominated by Si II and Fe II, 
however there are no spectra of 5pllAxl prior to day —31 
to conclude whether this is the case. 


3.2.5. Other SN Ib/Ic 


Plotted in Figure are the spectra of several SN Ib, 
SN Ic, and S N Ib/Ic^ke events, e.g ., SN 1984L, 1987M, 
and 2013ge ( |Drout et al. 2015). Similar to Fig- 


modeTspectra from 

Dessart et al. 

2012 

2015b 

I. 


in the middle panel ot h'igure [7 

we see the day —34 

Bmf5p09x4 model spectrum from 

Dessart et al. 

(2012) 


ward offset in wavelength for the 6250 A feature is ev¬ 
ident. Events that are similar to SN 2004gq include 
SN Ib 2004fe, 2009iz, 2009jf, and 2006ep during pre¬ 
maximum epochs. For the SN Ib 19901, the feature near¬ 
est to 6250 A on its red-most side is likely blended with 
a signature of He I A6678, and possibly with some con¬ 
tribution from C II. 

For the spectrum of the SN Ib 1984L, a He I A5876 fea¬ 
ture centered about ^ 5680 A implies a mean projected 
Doppler velocity of 10,000 ± 500 km s“^. Looking red- 
ward to the feature centered about 6250 A, associating 
the minimum with Si II A6355 implies a mean projected 
Doppler velocity of 5000 ± 500 km s“^. This differs from 
that inferred from He I by as much as 5000 km s“^. This 
suggests an association between Si II A6355 and the min¬ 
imum of the 6250 A feature is unlikely. 

Comparison between the day +8 spectrum of SN 1984L 
and the day —36 spectrum of the SN Ib 6p5Axl model 
from Dessart et al. (2015b) reveals moderate agreement. 
However, this particular benchmark model is not well- 
matched enough to invoke an unproven detection of 
Si II A6355. Specifically, the projected Doppler veloc¬ 
ities of the model and observations are not in sync. A 
visual comparison to the subsecment 6p5Axl model spec¬ 
tra (top-right panel of Figure 0 does not also reveal a 
more accurate match; i.e. the relative strength of the 
He I signatures improve, but the projected Doppler ve¬ 
locities remain too low. We emphasize that the 6p5Axl 
model has not been built for SN 1984L, yet this model 
still does not provide accurate predictions for any SN Ib 
shown in Figure This may be due, in part, to signa¬ 
tures of Si II A6355 that are too strong for the model. 

There does exist fair visual agreement between 
the 6170 A feature in the SN Ic 19 87M and the 
9000 km s“^ Si II profile produced by Dessart et al. 
(|2012|) for SN Ib 2008D; i.e., in our Figure [Tj overlap 
with the computed profile of Si II is arguably convincing. 
However, a composition that is contaminated by trace 
hydrogen is difficult to rule out since the epoch of the 
only available spectrum day +9) is not early enough 
to disaffirm tentative detecti ons of weak signatu res from 
hydrogen and helium atoms (Jeffery et alj|l991 ). 


For the S N Ic 2007gr, which is plotted in Fig¬ 
ures and [DessartjeF^ (2012) draw association be¬ 
tween the feature centered about 6350 A as a composite 
of Si II on the blue wing forming with a line velocity of 
^ 6000 km s“^ and Fe II near the minimum and red- 
ward. For the full fit, the day —32 Bmf5p09x3 model 
spectrum at these wavelengths is promising when com¬ 
pared to the day +13 spectrum of SN 2007gr. However, 
given the overall blue-ward displacement of the model 
compared to observations, we interpret the correspond¬ 
ing line velocities as being too high. 
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Rest Wavelength [A] 


Figure 7. Same as Figure 111 except now for select S N Ib/Ic. Spectrum ref erences: SN 1984L ||Wheeler & Levreault||1985|l; SN 1987M 
lIFilippenko et al. 1199 01 PefFe^ et al.ll991 ll: SN 19 901 llElmhamdi et a l.l2004ll: 200 4fe. 2004ga. SN :^U Ubge. blM ^UUYgr I telluric corrected), 
Sin 2UU91Z, HN 2UU9jt pJiodjaz et aI.|2U14| see also|Valenti et aI.|2Uil[|; SIN 2U13ge ||Drout et al.|20T^. 


Validation of this model for SN 2007gr could benefit, 
e.g., by a modest re duction of the kinetic energy released. 
(This was done by Branch et al. 2002 who used SYNOW 
to interpret and measure he II lines m the SN Ic 1990B.) 
However, the Si II would still remain too blue to confi¬ 
dently associate Si II with the absorption minimum of 
the 6350 A feature. If the computed spectrum were well- 
centered about this 6350 A feature, we suspect the ab¬ 
sorption minima between the observed and synthesized 
spectrum would be significantly offset at other wave¬ 
lengths. 

In the bottom of Figure one can also see that the 
computed spectrum produces a Si II A5972 feature that is 
not observed for SN 2007gr in this region, which is pos¬ 
sibly contaminated by narrow Na I. In addition, there 
is less of a “Si 11” notch on the blue wing of the larger 
6350 A feature to constrain the model when the ob served 
spectrum is co rrected for telluric features (cf. |Mod- 


jaz et al. 2014). Subsequently, we find the detection of 


Si II A6355 in SN 2007gr to be either unlikely and there¬ 
fore not representative of “Si II features” of SN Ib/Ic in 
general, or contribution from Si II A6355 is minimal for 
these kinds of core-collapse events. 


3.3. Super-luminous Supernovae 

The same distinctions used for classical types I and II 
are also used for super-luminous supernovae. SLSN II 
show conspicuous signatures of hydrogen, how ever the 
hydrogen fea ture is so far like that of SN Iln (Dessart 
et al. 2015aI, i.e. narr ow emission signatu res oi hydro¬ 
gen, el^, SN 2006gy (iSm ith et al.||2007|) and 2008am 
(Chatzopoulos et al. 2UII). 


For SLSN I, a spectral feature in the vicinity of 
6000 - 6400 A is either truly absent (or faint under 
noise), or a large and likely blended P Cygni profile is 
present. SLSN I spectra with a conspicuous 6250 A 
feature, i.e. similar to those observed for SN 2007bi, 
LSQI2dlf, and iPTFI3ehe, are frequently declared to be 
without a hydrogen feature, or even so far as “hydrogen- 


free” (Gal-Yam et al.|2009 Nicholl et al.|20I5 Yan et al. 


However, recall from Fig. I of Filippenko (1997) where 


the 6250 A feature in SN Ib/Ic spectra has been im¬ 
precisely labeled to signify an identification of “weak 
Si H” despite a conflicting definit ion of no silicon detected 
within §2.1 of Filippenko (1997). Under the assumption 
that “type I means no hydrogen,” an alternative inter¬ 
pretation of “weak Si II A6355” is often applied to the 
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Rest Wavelength [A] 

Figure 8. Model spectra from |Dessart et al.| (|2012|l (in orange) 
compared to the SN Ib 2008D and the SiN Ic 2UU7gr. i he blue Si II 
profiles computed for SN 2007gr and SN 2008D correspond to mean 
projected Doppler velocities of ~ 6000 km s“^ and 9000 km s~^, 
respectively. 


relatively large 6250 A P Cygni profile in the spectra of 
some SLSN I and all BL-Ic. 

Hence, if SN Ib/Ic and SLSN I shared the same “weak 
silicon” feature, then perhaps one would expect to see 
other lines of Si II in the spectra of BL-Ic and SLSN I 
where the feature is much stronger. Thus far, however, 
no strong evidence favoring a detection of other lines of 
silicon, e.g. Si II AA3858, 5972, has surfaced, while the 
low signal-to-noise ratios and low follow-up frequencies 
for these events prohibit more detailed investigations. 

In the following subsections we focus our attention on 
Figure where we have plotted comparisons between 
computed and observed spectra of BL-Ic and SLSN I 
found throughout the literature. Inconsistencies of mod¬ 
els that produce signatures of Si II A6355 are discussed 
in §3.3.1, and we again highlight spectral influence from 
Ha, or ejecta contaminated by unburned hydrogen, as a 
promising alternative in §3.3.2. 


3.3.1. Sill X6355 as a candidate identification 

Close inspection of Figure reveals an insufficient 
match between models that prMuce strong signatures of 
Si II A6355 and observations of BL-Ic and SLSN. A mis¬ 
match of “too blue” in wavelength is observed for a ma¬ 


jority of thes e fully-stripped C-1-0 models from Iwamoto 


et al. (19941, and either without any improvement in 


subsequent adaptations to BL-Ic and SLSN, or with the 
same conclusion of a model with a relatively shallower 
density profile compared to the observations. 

For the BL-Ic 1998bw the emission component pro¬ 
duced by the model near 6500 A is nearly 100 A too 
blue to match observations without considerable absorp¬ 
tion via ad hoc prescriptions of C II, Ne I, and Fe II. 


However, if the model were improved near 6200 A with 
C II and Ne I, these ions would be at odds with the dat a 
at other wavelengths (cf. Elmhamdi et al. 2006 2007). 
Alternatively, and using SYNAPPS, one could over-ht the 
region near 6300 A with C II and Fe II to secure a plau¬ 
sible match. Ultimately, however, the fit w ould suffer 
at blu er wavelengths, as shown recently by [Toy et ^ 
(2015) for SN 2013dx, which is similarly associated with 
a gamma-ray burst event. 


As for BL-Ic 2002ap ( Fole^_et_^|^ 
9000 km s“^ Si II profile from IDessart et a 

n 


2003), the 


(2012) re¬ 
veals similar problems of “too blue" without an ad hoc 


prescription of adjacent species that are not immedi¬ 
ately detected elsewhere. In particular, the discrepancy 
between the curvature of the computed profile domi¬ 
nated by Si II A6355 and that from observations indi¬ 
cates that Si II A6355 cannot be directly associated with 
the minimum of the 6250 A feature in BL-Ic 1998bw and 
2002ap. Similar conclusion s apply to so-cah ed “Hyper- 
novae”, such as SN 20 03hd (Deng et al. 2005) and 20031w 


(Mazzali et al. 


2006). In addition, the spectrum com- 


puted for SJN 20b3dh conflicts with observations through¬ 
out, including wavelength regions not shown in our Fig¬ 
ure [HI 

Comparison between the the March 7 spectrum of 
SN 2010ah (P TFlObzf) and the computed spectrum of 
Mazzali et al.|(2 013|) is shown at t he bo ttom of Figure 

Tc, i Mazzali et ^ (|2013 ) utilize rescal^ 
9981) to interpret spec- 
he 


For this BL-I 

C0138 models of|Iwamoto et al. ' 
tral signatures. Here again, we see the model does not 
provide adequate predictions. However, unlike previous 
comparisons to BL-Ic and SLSN, this is the first time in 
Figure]^ that we see a computed spectrum under-shoot 
the observed 6050 A feature. There are at least two av¬ 
enues for addressing the discrepancy between a profile 
dominated by Si II A6355 and observations of SN 2010ah. 


The blue wing of the 6150 A feature in the model 
conflicts with the observations, and similar to what 
is seen for SN 1997ef, 1998bw, 2002ap, 2003dh, and 
LSQ12dlf. One could posit that the model predicts 
velocities that are not as high as the observations, 
and that a simple scaling of kinetic energies, or a 
shallower ejecta density profile, would address the 
mismatch, thereby shifting the computed spectrum 
blue-ward. This is marginally true within the span 
of wavelengths plotted in Figure]^ However, such 
a correction for features near 6000 — 6400 A would 
create additional mismatch already seen for Ca II 
and Fe I I features at other w avelengths (shown in 
Fig 6. of Mazzali et al.||2013 ). 


- Alternatively, one could argue the mismatch be¬ 
tween the model spectrum does not correspond to 
the epoch needed for an appropriate comparison 
to SN 2010ah. However, had follow-up observa¬ 
tions been obtained, i.e. for more than the two 
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BL-Ic, GRB, SLSN-I/R 
Model Si II vs. Ha 



LSQ12dlf (SLSN-I) 
N14, B93 


SN2002ap (BL-Ic) 
D12 OSD, B93 


SN1998bw (BL-Ic) 
NOO, B93 


SN1997ef (BL-Ic) 
MOO, B93 


SN2007bi (SLSN-R) 


SN20031W (BL-Ic) 
M06 

SN2003dh (BL-Ic) 
M03, M06 

SN2010ah (BL-Ic) 
M13 


5600 6000 6400 6800 

Rest Wavelength [A] 


Figure 9. Comparisons bet ween compu ted spectra and a s election of B L-Ic, SLSN, SN 1 993J. and S N 2011fe. Spectrum references: 




(MUU) [hTTazzali < 


^_ .Idz UDUl): (MUij) rMazzali et aL|(|2003|l; fM06) iMazz ali et al.l 2UUbl) : fD12 UbU) b i II prohle tro r n a co mputed spectrum 

by |Dessart et aL] | [^ 1^, initially compared to tne blN lb 2008D; (iVllb) |Mazzali et al.] l |2013[ |; (N14) [Nicholl et al.| | |20r^ . 
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spectra obtained for SN 2010ah, we suspect the 
C0138 model spectrum would appear “too blue” 
near 6250 A at an ostensibly appropriate epoch. 

In Figu re the cornpute d spectrum for the BL- 
Ic 1997ef ( M^zali et al.|[M00 ) looks promising near the 
absorption minimum of the 6200 A feature. However, 
viewed as a whole, it seems the pa rent models of C021 


(Iwamoto et al. 19941 and C0138H (Iwamoto et al. 19981 
associated with the issortment of rnismatched synthetic 
spectra are better suited for other supernovae yet to be 
observed. 

The “type-R” SN 2007bi has been suggested to be pow¬ 
ered by radioacti ve decay through pair instability (Gal- 


Yam et al. 


2009). This origin , however, has since 


been 


contestecTl JVLoriyaetab 2CilO| Yoshida &: Umeda 2011 


Uessart et al. 2013 ' Ghatzopoulos et al. 2015). Since a 
spectrum has not been computed tor SiN 2007bi, one can 
only guess at its composition by comparing SN 2007bi 
to other objects in Figure Based on the similarity of 
6000—6400 A features, the 6250 A feature for SN 2007bi 
is either not dominated by Si II, or the evidence favoring 
a detection of Si II is insufficient. A similarly fleeting 
“identification” of He I (indicated by the arrows in Fig¬ 
ure]^ might instead favor either contamination from Ha 
near 6200 A, or a composition that is not necessarily 
hydrogen-free. 

Without sufficient improvements from stripped- 
envelope mode l s that produ ce signatures of Si II A6355, 


Nicholl et al 
to identity Si 


hi p 

ii ir 


2014, 2015) recently utilized SYNAPPS 
m the spectrum of another SLSN I, 
LSQ12dlf (shown near the top of our Figure |^. The 
interpretation of Si II as the 6250 A feature in LSQ12dlf 
is clearly inconsistent with the data from the red-most 
side of the emission component, through the absorption 
trough, and near the blue-most wing. Considering the 
observed feature slopes below the emission component of 
the proposed fit, which includes some Fe II, the mismatch 
is unavoidable without significant contribution from Fe II 
that, in the end, is difficult to constrain at other wave¬ 
lengths. M oreover, and similarly for SN 2013dx (Toy 
et al. 2015), the spectrum of LSQ12dlf is too noisy to 


assess the underlying structure of the composite feature 
near 6250 A. 


lYan et al.l (|2015[) recently presented observations of the 


to remain for massive stars, possibily as massive as 

- 125 M©? _ _ _ 

More recently, Arcavi et al. (2014 |2015 ) presented ob- 
servat ions of the rapidly-rising F'lF'lOiain. [Arcavi et al.| 
(2015) claim that the identification of its 6200 A feature 
can be either Si II A6355 or Ha. However, an interpreta¬ 
tion of “either-or” is not applicable to unproven identifi¬ 
cations of co mposite featur es. 


Curiously, Arcavi et al. (2015) find a fair match be¬ 
tween PTFlOiam and the 




I 1999em, apart from 


SN H-like H-Balmer lines that are absent for PTFlOiam. 


Arcavi et al. (2015) also compare the day -1-28 spec- 
truiu ot P'i'FlUiam to the day -1-11 spectrum of the 
SN la 1999ac, and claim that “the SN la fit is able 
to match this feature as Si H, as well as other fea¬ 
tures in the spectrum, with the major difference being 
that PTFlOiam has additional broad hydrogen emission 
lines.” 

However, we find another major difference is that the 
6200 A feature in PTFlOiam is not well-matched to the 
6135 A feature of SN 1999ac. In particular, the feature 
for SN 1999ac is noticeably too blue compared to that 
of PTFlOiam. Furthermore, if the SN la 1999ac can 
be used to identify the 6200 A feature of PTFlOiam as 
Si II A6355, while also “matching other features”, why 
then is PTFlOiam not classified as a SN la interacting 
with a H-rich mediumil^ 

The spectra of PTFlOiam are also noisy, and there¬ 
fore do not have a well-d efined minimum for the fea¬ 
ture near 6 200 A (see also Greiner et al.||2015 Leloudas 
et al.|[2015). However, for a 6200 A feature to be identi- 


hed as Si II A6355, this corresponds to a line velocity of 
8000 km s“^ (± 1000 km s“^). For the SN la 1999ac 


spectrum that Arcavi et al. (2015) claim to have matched 
to PTFiam’s 6200 A feature, its Si II A6355 line velocity 
is approximated to be ^ 11,000 km s“^. Consequently, it 
is unlikely that the 6200 A feature of PTFlOiam can be 
solely attributed to Si II A6355, and we suspect models 
that produce spectra dominated by Si II A6355 will be 
unable to unequivocally match both the minimum and 
the blue and red wings of the 6200 A feature without 
producing inconsistencies elsewhere. 

This reliance on Si II A6355 as a candidate identifi¬ 
cation reveals an ongoing problem. So far, there has 
been no observational preced ent for lone signatur e s of 


SLSiN I iP'iF'13efie that rev ealed strong emiss ion in Ha Si II A6355 outside of SN la ([Wheeler et al. 


this CSM material originates from the progenitor in the 
years prior to the explosion, and estimate a CSM mass 
of < 30 M©. Given the progenitor masses for some of 
these SLSN I are thought to be 67—220 M©, this would 
imply the progenitor’s mass is relatively high, and pos¬ 
sibly with trace amounts of hydrogen by the time of the 
explosion. (This of course depends on the episodes of 
mass loss throughout the progenitor’s evolution.) 
Without providing a compute d spectrum that is suf- 


1995 Fil- 


during late-nebular phases. Yan et al. (2015) conclude ippenko 1997), nor has an accurate search ot additional 


ficiently matched to iPTF13ehe, Yan et al. (2015) also 
interpret the 6000—6400 A feature as predominately 
Si II A6355. Given the evidence for interaction between 
the ejecta of this SLSN I with material contamin ated by 
hydrogen that was previously ejected (see also Chugai 


& Chevalier 2006 Milisavljevic et al. 2015), should we 


then expect trace amounts of hydrogen (~ 10 ^ M©) 


signatures of Si II been carried out. Yet the common 
wisdom for SN Ib, Ic, BL-Ic, and SLSN continues to de¬ 
fault to a subjective and inconsistent identification of 
Si II A6355. It is therefore a wonder why Si II A6355 has 
be en so strongly fav ored as an int erpretation (cf. Fig. 1 
of Filippenko 1997 and Fig. 2 of Modjaz et al. 2014). 
A more likely alternative identification is not necessarily 
Ha, much less Ha alone. Rather, a 6200 A feature that 
is influenced by trace amounts of ejected hydrogen would 
be supported by the same observations of PTFlOiam pur- 

3 PTFlOiam is not likely a SN la because it does not look like a 
SN la prior to day +28. Yet an identification for its 6200 A feature 
is being interpreted as Si II A6355 based on an unmatched com¬ 
parison to the SN la 1999ac, while an identi fication of Si II A6355 
is so far only de finitely detected for SN la dWheeler et al.[[1995[ 
[Filippenkolig^. -- 
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Days Since Maximum Light 


Figure 10. Estimated line velocities of SN Ib 2008D from its 
absorption minima. Vertical error bars shown for Si II and Fe II 
represent a minimum error of it 750 km s“^. The errors could be 
as high as 2000 km s“^. Sh own above are values obtained for the 
SN la 201 Ife using SYNAPPS ([Parrent et al.|2012[|. 


portedly in teractiiiK with a H -rich medium near maxi¬ 
mum light (Arcavi et al.||2015|. 


3.3.2. HI \6563 as a candidate identification 

Compared to Ha in SN II spectra, SN Ilb show weaker, 
albeit typical, sign atures of hydroge n alongside conspic¬ 
uous helium lines (Filippenko 19881. The Ha profile for 


SN Hb is also generally broader than that for SN H be¬ 
cause the mass of the remaining hydrogen is greater for 
SN H than SN Hb. Consequently, the relative weakness 


of Ha in SN Hb spectra enables adjacent blending from 
species such as Si II, Fe II, and possibly higher velocity 
H I to broaden a spectral feature that would otherw ise 
be solely attributed to Ha (Milisavljevic et al.||2013 and 
references therein). 

Considering the similarity of all BL-Ic and select 
SLSN I/R plotted in Figure we find contamination 
by Ha to be a plausible solution. One caveat is that 
the classically defined rest-wavelength P Cygni emission 
component of the observed feature is no t centered about 
a rest- wavelength of 6563 A. For SN H, Anderson et al. 
(20141 have shown the blue-shifted emission of Ha is sen¬ 


sitive to the ejecta density profile, particularly when the 
radial falloff is steep. Blue-shifted emission components 
in Ha have also bee n observed for SN H b, e.g., the well- 
observed SN 1993J (Baron et al. 1993|). Therefore, an¬ 
other solution that would support an interpretation of 
hydrogen for BL-Ic and select SLSN I/R is if the un¬ 
derlying emission component of Ha is blue-shifted along 
with the corresponding absorption minimum. 

3.4. Measures of projected Doppler velocity 

Estimates of so-called detachment velocities for hydro¬ 
gen are either accurate from first principles, or direct 
measurements are again spoiled by effects of radiation 
transport in that the absorption minimum does not ac¬ 
curately estimMeji^n projected Doppler velocities for 
Ha (cf. Branch |1977). Under the assumption of either 
localized line formation, or multiple resonance line scat¬ 
tering (effectively SYNAPPS), line velocities inferred from 
absorption minima may over-shoot the optimum posi¬ 
tion of these ostensibly found hydrogen-deficient regions 
of ejecta. If faint Ha in SN Ib, Ic, BL-Ic, and SLSN I/R 
behaves similarly compared to SN H and Hb, then trace 
amounts of detectable hydrogen would indicate high pro¬ 
jected Doppler velocities in spite of being relati vely closer 


to the pho tospheric line forming region (cf. Anderson 

eFallpnt . 


By contrast, the red-shift needed to rectify ill-matched 
signatures of Si II A6355 is on the order of 200 A. The 
effect that opposing interpretations of Si II A6355 ver¬ 
sus H g ha s on estimates of line velocities is shown in 
Figure [T0| for the SN Ib 2008D. Assuming one line domi¬ 
nates the 6250 A feature, the error would be on the order 
of 5v ~ 9000 km s“^ between Si II A6355 and Ha; i.e., no¬ 
ticeably large when spectroscopic resolutions of blended 
features are on the order of > 750 km s“^ for these classes 
of SN 1. Compared to the SN la 201 Ife, where signatures 
of both Si II and Fe II can be used to trace respective 
line velocities, the same conclusion is not reached when 
considering the SN Ib 2008D. 

Similarly, line velocities of Fe II are also noticeably 
higher than that of Si II A6355 for SN 19941, 1998bw, 
2002ap, 2004aw, 200 9bb, and 2012ap. show n in Figure [Tl^ 


(see also Table 2 of Stathakis et al. 
continues into post-niaximum epochs 


2000). This trend 
hat are not shown, 


and indicates that if a signature of Si II A6355 is shaping 
the spectrum, then it is weak and contributing to the 
blue-most wing of 6000 A absorption features (instead of 
dominating an entire profile, much less the minimum). 

4. SUMMARY AND CONCLUSIONS 

In this work we examined publicly unavail¬ 
able model spectra from both “hydrogen-poor” 
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Days Since Maximum Light 


Figure 11. Estimated line velocities of 6000 A features in SLSN 
and BL-Ic spectra. The grey points denote values obtained un¬ 
der an assumption that the absorption minimum is identified with 
Si II A6355. Estimates of line velocities for Fe II are the points 
shown in black. For some of the available observations, we ex¬ 
cluded points that did not have either an obvious minimum or 
underlying structure. 

(Mjj ^ 10“^ — 10“^ Mq) and “hydrogen-depleted” 
(Mjj ~ 10“® Mq) compositions and compared them to 
observations of multiple type I supernova subclasses. 
Apart from natural uncertainties from line blending, 
and regardless of the spectrum synthesizer used, all 
1-dimensional model spectra intended for thermonuclear 
SN la are well-matched (Figure [^. Since the evidence 
favoring a detection of Si II A6S55 is overwhelming, 
direct association between Si II A6355 and the minimum 
of 6150 A features can be made in practice. 

However for SN Ib, Ic, BL-Ic, and SLSN I/R, the evi¬ 
dence favoring detections of Si II A6355 spans few events. 
This implies the direct association between Si II A6355 
and the minimum of 6000—6400 A feature is not accu¬ 
rate. In addition, the observational evidence that would 
otherwise support an interpretation of Si II A6355, e.g., 
signatures from lines other than A6355, is as fleeting as 


any weak evidence that would similarly confirm the pres¬ 
ence of trace hydrogen through H/I (Figure 

Among the available spectra stemming from models 
of stripped-envelope supernovae, and compared to ob¬ 
served 6000—6400 spectral features of SN Ib, Ic, BL-Ic, 
and SLSN, we find a consistent blue-ward mismatch in 
wavelength, in part, due to signatures of Si II A6355 that 
are too strong. Subsequent explanations for this phe¬ 
nomenon have included a shallow density pro file for the 
corresponding models (|Hachinger et al. 2012| and refer¬ 
ences therein), or a prescription that is out of sync with 
the time-series observations (see our §3.3.1). However, 
these notions are counter to both the successful interpre¬ 
tation of Si II for SN la (Figure]^, and the categorical 
shortcomings of stripped-envelope models that produce 
signatures of Si II that are too strong for pre-maximum 
observations of SN Ib, Ic, BL-Ic, and SLSN. 

The inability of stripped-envelope models to accurately 
reproduce the shape of a given spectral feature could 
be attributed to many potential issues with velocities, 
composition, excitation/ionization, and assumptions of 
spherical symmetry in the model, and would not neces¬ 
sarily imply that a given model is fundamentally wrong 
(or right). However, if both photometry and time-series 
spectra of a given event are to be thought of as sensitive 
to the radiation field and properties of the ejecta, then 
a model that is only approximately consistent with one 
of these pieces of information precludes firm conclusions 
about the nature of the spectral sequence being studied. 

Contribution from species with line signatures near 
6250 A, such as C II and Fe II, may help to both pre¬ 
vent a poor match for a model a nd refine prescription s 
for empirical measurements (cf. Dessart et al. 2015b). 
This solution, however, will require simultaneous agree¬ 
ment between the observed and computed photometry 
since this is so far not the case when an interpretation of 
Si II with Fe II has been put forward (see §3.2). 

Furthermore, today’s spectrum-limited observations 
do not provide the necessary constraints for invoking de¬ 
tections of C II through both A6580 and A7234 coun¬ 
terparts. This underscores the need for additional sur¬ 
veys that embark on maximizing follow-up frequency 
and signal-to-noise ratios throughout the evolution of the 
spectrum. 

Among the most promising matches with stripped- 
envelope core-collapse models, two “realistic” model 
spectra can be described as being in fair agreement with 
two observations from two different objects; namely the 
hydrogen-poor (M// ^ 10“^ Mf?^) and compositio nally 
mixed SN Ib model Bfm4p4Ix4 (Dessart et al. 2012) com¬ 
pared to the SN Ib 1999dn, and the hydrogen- depleted 


(Mff ^ 10 ° M 0 ) SN Ib model Bmf5p09x4 (Dessart 


et ai.j2012 ) compared to the SN Ib 2004gq. Even so, the 


rise-times for these models are not in good agreement 
with observations fo r these and most other SN Ib/Ic (cf. 
Wheeler et al.|[2M5), which brings into question the va- 


lidity of the models altogether. 

One might expect detailed modeling through PHOENIX 


^ Although section 2.1 of|Filippenko|(|l997[ claims Si II A6355 is 
detected for SN la, but not tor blN Ib/lc, th e labels of Si II A6 355 
against SN la, Ib, and Ic spectra in Fig. 1 of |Filippenko| (|1997|| are 
at odds with this statement. As a result, reterences made to this 
review have sided with either the labels of Fig. 1 or the notions 
expressed within the body of the text. 
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and CMFGEN, both of which assume spherical symme¬ 
try for post-processing, to show some improvement at 
wavelengths red-ward of 6250 A for hydrogen-depleted 
compositions since these synthesizers at least consider 
a larger volume of line formation than those methods 
that assume a sharp photosphere. However, such a red- 
ward fix does not appear to reveal itself in either the 
PHOENIX or CMFGEN synthetic spectra displayed in our 
Figures and[^ In addition, an apparent absence 

of 200 AHine-shifts toward redder wavelengths for other 
spectral features appears to suggest that considerations 
of 3-dimensional radiation transport would not conspire 
to enable the minimum of a Si II A6355 profile to coincide 
with the observed minimum near 6250 A. 

Thus, an absence of faint Ha originating from trace 
hydrogen (~ 10“^ Mq) is not simply a given for these 
“type I” supernovae, while the common wisdom has been 
to judge the quality of a computed spectrum from a fully- 
stripped model by its [in] ability to match an unproven 
signature of Si II A6355 within the feature centered about 


6250 A. Moreover, observations of would-be conspicu¬ 
ous signatures of hydrogen in the spectra of BL-Ic and 
SLSN I suggest that the broad emission component of 
Ha undergoes a significant blue-ward shift relative to the 
rest wavelength of 6563 A, which is an effect known to 
arise from steep dens ity-gradients in the outermost ejecta 


of SN H in general (Baron et al. 1995 Anderson et al. 


20I4I. 

In the literature there are at least three promising 
models with computed spectra where ^ 10“^ Mq of hy¬ 
drogen is present within the outermost layers of ejecta, 
and the times-series spectra for only one of these mod- 


(Dessart et al.|20I5 

bl), the aforementioned Ha-negative 

SN Ib Bmf4plix4 

Dessart et al. z 

012), 

and the cal- 

culations done by 

ames & Baron 

m 

h-f 

o 

) who previ- 


ously showed that Ha cannot be ruled out assuming 
extended regions of trace amounts of hydrogen for the 
SN Ib I999dn. Hence, whether or not the computed spec¬ 
trum for a hydrogen-poor event is Ha-positive, the mass 
of hydrogen in the outermost layers of some caught-early 
SN Ib, and therefore the progenitor, need not fall be¬ 
low ~ 10“^ Mq for the model to remain consistent with 
6250 A spectral features (§3.2.1). 

Depending on the sample of objects utilized to infer the 
progenitor and mass-loss history of “hydrogen-poor” su¬ 
pernovae (e.g., Hb/Ib), relationships may exist between 
properties of the light curves and the extent of trace hy¬ 
drogen throughout the ejecta. Given a suite of improved 
model spectra, faint detections of Ha, as well as non¬ 
detections stemming from trace hydrogen in the ejecta, 
might also be used to gain insights into the distribution 
of hydrogen near the surface of the progenitor star, or 
binary product, at the time of explosion. 
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JTP is indebted to David Branch and J. Craig Wheeler 
for extensive discussions and valuable advice during the 
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Table 1 

Common Lines in Optical Spectra of Type I Supernovae 


Ion 

Rest Wavelength (A) 

“C 11'’'= 

6580, 7234 


7774, 8446, 9264 

“Mg 11*’“ 

4481, 7896 

“Ca 11*’“ 

3969, 3934, 8498, 8542, 8662 

“Si II 

3838, 4130, 5051, 5972, 6355* 

“Si III 

4560, 5743 

“S II 

4163, 5032, 5208, 5468, 5654, 6305 

“Fe 11*’“ 

4025, 4549, 4924, 5018, 5169 

“Fe III 

4420, 5156, 6000 


If the ‘a’ is on the left, an ion has been de¬ 
tected for SN la; “be” on the right, an ion has 
been detected for SN Ib and Ic. Fe III may or 
may not be detectable, if present, on account 
of severe line-blending for all SN Ib and Ic. 
For fast-evolving events, such as SN 2005hk 
and 2010X, an identification of Si II A6355 also 
seems likely (Drout et al. 2013). 

* When Si II is identified in the spectra of 
SN Ib, Ic, BL-Ic, and SLSN, it is based on 
a single line, A6355, i.e. unconfirmed in the 
same sense as for Ha. 



